Rapidity dependence of large-p^ hadron production at RHIC 
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We study the dependence of parton energy loss (quenching) on rapidity in ultra-relativistic nuclear 
collisions at RHIC. This can provides invaluable information on the density of the medium, which 
should be more dilute going away from mid-rapidity, thereby reducing the effect of quenching. We 
predict a clear effect at moderate transverse momenta ~ 3 GeV. 
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Energy loss of high energy quark and gluon jets in the 
dense matter produced in ultra-relativistic heavy ion col- 
lisions is expected to lead to jet quenching and can po- 
tentially probe the quark-gluon plasma formed in these 
reactions M |, |, | g |, |{. Recent preliminary data 
from the first run of RHIC at BNL suggest the existence 
of this phenomenon || , where moderate pr ~ 2 — 6 GeV 
hadron spectra were found to be suppressed by a factor 
^2 — 5 relative to the pQCD expectations. In addition, 
the suddenly saturated azimuthal asymmetric moment 
V2(pt) for values of pr > 2 GeV was interpreted as ev- 
idence of jet quenching in a dense gluon plasma with 
rapidity density of the initially produced gluons in the 
range of dN^/dy ~ 500 - 1000 [§. However, to finally 
constrain the amount of jet quenching which affects high 
momentum hadrons at RHIC, we must seek further in- 
dependent observables to test these phenomena and sup- 
port their existence. 

Previous studies mostly focus on the suppression of 
Pt spectra in nucleus-nucleus collisions compared with 
the proton-proton case at the same energy, or in central 
nucleus-nucleus collisions compared with the peripheral 
ones. This is always done at mid-rapidity. The point we 
want to stress is that energy loss effects depend on the 
density of the environment, i.e., higher medium densities 
mean higher amount of energy loss. Although one can 
study the dependence on the medium density by vary- 
ing the impact parameter, at fixed impact parameter b 
for nucleus- nucleus collisions, also different rapidities pro- 
vide different densities of the medium, through which the 
high momentum jet travels. As illustrated in Fig. 1, when 
a high energy jet travels through the quark-gluon plasma 
produced at the early stage of nucleus-nucleus collisions, 
the medium density is much larger at mid-rapidity than 
in the forward or backward directions (large rapidity). 
Because jet energy loss is proportional to the density of 
the local medium^, [| 0, ||, |[ |?], ^ , at different rapidities 
the energy loss of the fixed pr jet will be different, and 
the large pr hadron spectrum from jet fragmentation will 
also have different behavior. Therefore, by scanning the 
rapidity of hadrons produced in central nucleus-nucleus 
collisions, one can observe the variation of energy loss 
effects and, as we will show with the following calcula- 
tions, its influence on the rapidity spectrum of charged 
hadrons at fixed pt- A study based on an analogous idea 
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FIG. 1: Illustration of the spatial extention a nucleus-nucleus 
collision at fixed time in the cm. frame. Dashed lines show 
trajectories of jets produced with small (a) and large (b) val- 
ues of y, traversing dense and dilute parts of the medium, 
respectively. 



was recently made for J/ip suppression as function of xp 
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Large px hadron production in pp collisions can be 
expressed in a standard pQCD form as 

da h f 

= K I dx a dx b d 2 k a d 2 k b g p (k a ) g p (k b ) (1) 



dyd 2 p T 

X fa/p{Xa,Q Z ) fb/p(xb,Q Z ) 



1 , r - da a b^cX D h / c {z c ,Q 2 ) 

di 7TZ C 



where f i i p is the distribution function of parton i in the 
proton, da/dt are parton-parton coss sections and -D/w c is 
the fragmentation function of parton c in hadron h. We 
include intrinsic transverse momentum effects adopting 
the phenomenological /cy-kick model with the Gaussian 
distribution for the intrinsic kx of a parton inside the 
proton 



g p (k T ) = (ttA^)- 1 exp(-£4/A^). 



(2) 



In the numerical calculations, we choose the GRV94 pa- 
rameterization for the parton distribution functions jll| , 
and the KKP set for the fragmentation functions into 
charged hadrons With the scale choice of Q 2 = p 2 , 
and an overall factor K = 1.5, we can reproduce charged 
hadron spectra in pp collisions in various experiments 
and also the UA1 data for pp collisions at y/s = 200 



2 



GeV. There we use the energy dependent values for A 2 
H (A 2 = 1.0 GeV 2 at */~s = 200 GeV). 

We now turn to the case of proton- nucleus collisions. 
The first feature we need to discuss is that of parton 
shadowing. Although some approaches attempt to solve 
the problem in general terms, a complete solution is not 
yet at hand. Since this aspect is not our main focus, 
we make use of the popular EKS98 parameterization Q 
and modify it to take into account impact parameter de- 
pendence. This leads to the factor 



R i/A (x,Q 2 ,b) = 1 



T A {b) 
(Ta) 



S % / A {x,Q 2 ) ■ 



(3) 



where (T A ) = A~ x f d 2 b T\(b) is the mean value of the 
thickness function and Si/ A (x, Q 2 ) is the shadowing fac- 
tor of the EKS parameterization. In this way we satisfy 
the relation 

(R i/A )(x,Q 2 ) = A- 1 J d 2 bT A {b) R t/A ( Xl Q 2 ,b) 

= l-S l/A (x,Q 2 ) (4) 

which reproduces, by construction, the EKS98 parton 
distribution ratios. 

Another feature which appears in pA collisions is the 
broadening of transverse momentum spectra (Cronin ef- 
fect). This is understood as parton rescattering during 
the traversal of the projectile through the nucleus. This 
effect is parameterized by convoluting the spectra with 
a Gaussian function whose width is proportional to the 
mean number of rescatterings the projectile suffers inside 
the nucleus |H| . With this method, the pr broaden- 
ing effects on hadron production in pA and AB collisions 
have been previously studied with a modification of the 
intrinsic fcy distribution of a nucleon by convoluting it 
with the broadening factor Therefore, the ini- 

tial intrinsic kx distribution of partons will be modified 
with respect to Eq. (|2|) to the form 



g%(k T ,b) = (nA A (b))- 1 eM- 



-4/a»), 



(5) 



with 



A 2 A (b) = A 2 p + A 2 n (b), A 2 n (b) = Xa NN T A (b), (6) 

where A is a coefficient corresponding to the average value 
of the additional squared momentum, acquired in each 
collision by parton i, and aN^T A (b) is the number of 
scattering the projectile nucleon suffers inside the nu- 
cleus. In our calculations, we also take into account the 
broadening effects to the final jet spectrum prior to frag- 
mentation and convolute it with the broadening distri- 
bution 

g A m (k T ,b) = ( 7 rA 2 m (&))- 1 exp(-4/A 2 m (6)). (7) 

with A 2 in (b) = A 2 n (b). In the following we differentiate 
the broadening effects between quarks and gluons, where 



the quark one is 4/9 smaller than the gluon one. By fit- 
ting the charged hadron spectrum from various pA data, 
we obtained the coefficient A = 0.1 GeV 2 , which roughly 
agrees with previous studies |15, 
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Considering now the AB case, we first of all generalize 
the nuclear effects as discussed above. This is straight- 
forward for shadowing and initial broadening. In fact 
each colliding nucleus is affected and not only the target 
as in the pA case. Concerning final state broadening we 
must now take into account both nuclei at once, there- 
fore A 2 in = \<jnn{T a +Tb). Then we come to discuss 
the key issue of the present paper, that is parton energy 
loss. In the literature, there have been intense studies on 
this subject J|, [| [|, [|, ||, These studies show that the 
energy loss of the fast parton in a dense medium depends 
on the medium density, the initial parton energy and the 
path length of the parton traversing the medium. In dif- 
ferent kinematical regions, the total energy loss has a 
different functional dependence on the above mentioned 
factors |l7| . However, in the region of interest for the pro- 
duction of few GeV hadrons, it is widely believed that the 
energy loss is proportional to the parton density and to 
the square of the path lengthen a static medium). Then, 
we make a simple observation. If a jet is produced at mid- 
rapidity, it will encounter a dense medium, while at large 
enough rapidities the energy loss effect disappears. We 
therefore assume that a parton c loses energy because of 
gluon radiation induced by the medium, in the amount 



AE c (y,E c ) = q c E c 



dn 
dy 



(8) 



where q c is the average value of energy loss which implic- 
itly also includes the size of the medium. Its value has to 
be determined and for quarks it is 4/9 that of gluons. In 
the simple parameterization above, we assume that the 
energy loss is also proportional to the original jet energy, 
according to the recent findings within the GLV formal- 
ism H Q. We note that in the literature, there exist 
various formalisms which give rise to a different energy 
dependence of the parton energy loss || ^, |l7j. In the 
GLV formalism, the linear energy dependence of jet en- 
ergy loss adopted here is mainly caused by the kinematic 
limit of the gluon radiation at low energy Q . This lin- 
ear dependence may not be valid at large energy. In fact, 
at asymptotically high energies, the dependence becomes 
logarithmic (|, ||. Moreover, a recent study on this sub- 
ject shows that thermal gluon absorption can reduce 
the energy loss for a low energy jet and change the energy 
dependence of the energy loss to stronger than linear at 
low energy. Nevertheless, in our calculations we concen- 
trate on the production of few GeV hadrons where the 
linear dependence for jet energy loss should be a reason- 
able approximation ||, 0, |T^] . 

Eq. (||) also contains the rapidity dependence of jet 
energy loss, which relies on the density of the medium 
dn/dy. We note that at large rapidities the time evolu- 
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tion of the medium may be quite different from that at 
mid-rapidity (y = 0), and jet energy loss might have a 
complicated dependence on geometry. However, we ar- 
gue that at larger rapidity, the dominant effects to the jet 
energy loss is caused by the lower medium density. This 
means less energy loss at larger rapidities. So, as a first 
estimate of the rapidity dependence, in our calculations 
we made the simple but intuitive assumption contained 
in Eq. (||) to parameterize the rapidity dependence of jet 
energy loss. For the medium density dn/dy, we further 
assume for simplicity that the plasma density at the early 
stage is proportional to the measured rapidity density of 
charged hadrons, therefore including both soft and hard 
produced particles. The shape was fitted as in Jl(| to the 
recent PHOBOS data and here normalized to 1 at y = 0. 
We do not take into explicit account the time evolution 
of the plasma which, due to geometric and dynamical ef- 
fects, requires a more detailed study of its properties. We 
therefore consider only the pre-factor q c which contains 
the average effect of the space-time evolution. 

If a parton loses energy, its fragmentation is then mod- 
ified 1 19, since its momentum fraction which is con- 
verted into a hadron is reduced. This can be approxi- 
mately realized by introducing the momentum fraction 
zL= p±/ {p c - AE c (p c )), being p c = Pa_/ z c (Although see 
pffl), and changing the fragmentation function as 



D 



h/c 



(z c ,Q 2 ) 



D c / h (z' c ,i 



(9) 



In this way we neglect a small contribution due to the 
fragmentation of radiated soft gluons. 

Substituting all the elements which appear in describ- 
ing nucleus-nucleus collisions in addition to the pp case, 
we obtain the final form for the hadron spectrum 



dN h AB 
dyd 2 p T 



(b) = K J d 2 b A dx a dx b d 2 k a d 2 k b d 2 k c (10) 
x g l ^(k a ,b A ) R a/A {x 

x g B (k b ,b B ) R b/B (x b ,Q 2 ,b B ) h/ p (x b ,Q 2 ) 



da ab ^ 



cX JmiZ 



dt 



g J AB (k c -p T /z c ,b A ,b B ) 



z' c D h/c (z' c ,Q}) 



In this way we have completed the discussion of the tech- 
nical part of the paper and proceed to examine the results 
of the computations and their physical implications. 

Our numerical results are presented in Figs. 2-3, where 
we display ratios of spectra in nucleus-nucleus collisions 
respect to pp collisions. Actually we computed 



dN ^(b) 



J? h f„ r, h\ — d V d2 PT 
RaB\.ViPT,0) - 



dal 



(11) 



PT 



where T AB (b) is the usual overlap function, and took 
6 = 0. Fig. 2 is for the ratios of pr spectra at different 
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FIG. 2: Transverse momentum dependence of the ratio of 
spectra in AB collisions with respect to the pp case, for dif- 
ferent rapidity values. 
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FIG. 3: Rapidity dependence of the ratio of spectra in AB 
collisions with respect to the pp case, for different transverse 
momentum values. 



rapidities, and Fig. 3 for the ratios of rapidity spectra at 
different pr values. All the results are for Au + Au colli- 
sions at y/s — 200 GeV, relevant to RHIC. For the energy 
loss coefficient we use q c = 0.27 with which we can repro- 
duce the suppression of the charged hadron spectrum in 
central collisions observed at y/s — 130 GeV. In Fig. 2, 
we see that for different rapidities the ratios have differ- 
ent behaviors. In general, at larger y, the suppression is 
reduced, as expected, due to the decrease in the medium 
density. Especially, at very large rapidity one can re- 
ally observe the Cronin effect. The rise at pr > 7 GeV 
of the curve for y — 3 arises due to the vicinity of the 
kinematical boundary for hadron production, where our 
treatment of final state px broadening as folding with 
Gaussians is no more valid. 

Even more interesting results can be found for the ra- 
tios as functions of rapidity at different values of pr, as 
illustrated in Fig. 3. At mid-rapidity, because of the high 
density of the medium, the largest energy loss produces 
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the strongest suppression. As rapidity increases, the en- 
ergy loss effect dies out and the ratio increases. Espe- 
cially for the largest values of p? what is clear is the 
presence of a concavity around mid-rapidity in the ra- 
tio. This feature should not be present in absence of jet 
quenching, while our calculations provide it clearly. This 
is the result of our investigation and we argue that it 
should be a good signal to observe energy loss at work. 

The idea of scanning the medium by computing rapid- 
ity spectra can be extended to the study of the y depen- 
dence of the V2 coefficient of azimuthal asymmetry. We 
expect its value to vary significantly with rapidity. More- 
over, a comprehensive study of direct photon production 
could be a cross check of the present work, since pho- 
ton should not experience energy loss end therefore their 
spectra should not be quenched. 

In conclusion, in this paper we have calculated the ra- 
pidity dependence of hadron spectra for high values of 
Pt in Au + Au collisions at ^fs = 200 GeV. Compar- 
ing the results to the pp case, we found a strong signal 
to test the effect of jet energy loss. The observation of 
such a behavior within future RHIC data will give a new 
test and support for the idea of jet quenching, therefore 
providing important information about the existence and 
properties of the quark-gluon plasma. 
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